Mon. Not. R. Astron. Soc. 000,[l]|71(2004) Printed 2 February 2008 (MN MbX style file v2.2) 



The clustering of Luminous Red Galaxies around Mg ii absorbers 
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ABSTRACT 

We study the cross-correlation between 212 Mgll quasar absorption systems and ^ 20,000 
Luminous Red Galaxies (LRGs) selected from the Sloan Digital Sky Survey Data Release 1 
in the redshift range 0.4 ^ z ^ 0.8. The Mg II systems were selected to have AA2796 & 2803 
rest-frame equivalent widths ^ 1.0 A and identifications confirmed by the Fe II A2600 or Mg I 
A2852 lines. Over comoving scales 0.05-13/i~^ Mpc, the Mg II-LRG cross-correlation has an 
amplitude 0.67 ± 0.09 times that of the LRG-LRG auto-correlation. Since LRGs have halo- 
masses greater than 3.5 x 10^^ M© for Mr < —21, this relative amplitude implies that the 
absorber host-galaxies have halo-masses greater than ~ 2-8 x 10^^ M0. For lO^'^ M0 LRGs, 
the absorber host-galaxies have halo-masses 0.5-2.5 x 10^^ M0. Our results appear consistent 
with those of Steidel et al. ( 1994) who found that Mgll absorbers with W^'^sH ^ q.SA are 
associated with ^0.7L|j galaxies. 
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1 INTRODUCTION 



The con nection between quasar (Q SO) absorption line systems and 
galaxies dBereeron & Boissell99lh is important to our understand- 
ing of galaxy evolution. Absorption lines provide detailed informa- 
tion about the physical conditions and kinematics of galaxies out 
. to large impact parameters {R > 100 kpc), re gardless o f the ab- 
' sorbe r's intrinsic luminosity (e.g. iRauch et alj |1996: Ell ison et alj 
l200(i) . Mg II AA2796 & 2803 are amoungst the most studied metal 
lines since the doublet signature makes for easy detection. 

Past results show that Mg II absorbers are not unbiased trac- 
ers of galaxies but are biased towards late-type galaxies which 
do not evolve strongly from z ~ 1. The morpho l ogical con- 
straints come from imaging by ' Steidel & Sargend Jl992h and 
[steidel. Dickinson & Persson (1994) who found that Mgll ab- 
sorber host-galaxies have if-band luminosities consistent with nor- 
mal 0.7Lr Sb galaxies. Further Hubble Space Telescope imaging 
jSteideil993) indicated that Mg II absorbers at Zabs < 1 are drawn 
from field galaxies of all disk morphological types. These galaxies 
are also found to have roughly constant star formation rate since 
2 ~ 1: from a sample of 58 Mg II absorbers with re st-frame equiv- 
alent widths ■py.'^'sii ^ 0.30 A. ISteidel et alj <1994l) found that the 
mean rest-frame Mb and B — K colours of the host-galaxies do 
not evolve in the redshift range 0.2 < Zabs < 1- Furthermore, their 
rest-frame K-hanA luminosity function (LF) closely matches the 
/f-band LF at z = down to O.OSI/^- This 'no evolution' picture 
of the A'-band LF implies little or no stellar mass evolution from 
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z ~ 1. However, using more di rect constraints from deep imaging 
survevs- loTckinson et al]<2003h and lRudnick et alJ<200 3') find that 
the stellar mass increased by a factor of 2 over this epoch. Thus, 
these earlier studies imply that Mgll absorption-selected galax- 
ies are biased towards non-evolving, luminous disk morphological 
types. 

From the observed absorber-host-galaxy impact parameter 
distribution, Steidel ( 1993, 1995) constrained the cross-section of 
Mgll absorbers with iy,^'sn j> o.30A to have radius i?x ~ 
40/i~^ kpc (physical, 70/i~^ kpc comoving). In addition, these 
systems are always found to be associated with neutral hydrogen 
absorbers in the Lyman limit regime (log N-^ j 3 x lO^'^ cm~^). 

Little information currently exists about the environment 
of Mgll absorbers on scales up to lO/i^^ Mpc. Recently, 
iHaines. Campusano &. Clowes! d2004ft analysed the clustering of 
early-type galaxies around two Mg II absorbers at «abs = 0.8 & 1.2 
using wide field images (40' x 35')- They find a significant excess 
of galaxies across the field and conclude that large-scale structures 
containing Mg II absorbers mark out volumes of enhanced galaxy 
density. 

All the above studies are based on relatively s mall (< 50) 
samples of Mg II absorbers and, with the exception of lHaines et alJ 
(2004), relatively small-area (< 10' x 10') galaxy sur veys. The 
Sloan Digital Sky Survey (SDSS: IStoughton et aj2002h allows us 
to significantly transcend these limitations. In this paper we use 
~ 20,000 Lu minous Red Gala xies (LRGs ) from SPSS D ata Re- 
lease 1 (DRl: lAbazaiian et alj l2003: Stra uss et alj|2002h to con- 
strain the environment, and more specifically, the mass of the ha- 
los associated with 212 Mg II absorbers. In a hierarchical galaxy 
formation scenario, the amplitude ratio of the Mg II-LRG cross- 
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Figure 1. Example Mgll absorption system showing the Mgll 
AA2796/2803 doublet and the supporting transitions, Fe II A2600 and Mg I 
A2852. Note also the Fe II A2587 line at Aots ~ 4415A. The upper solid 
line is our fitted continuum and the lower solid line is the 1 cr eiTor array. 



correlation and LRG-LRG auto-correlation is a measure of the 
relative masses of the halos associated with Mgll absorbers and 
LRGs. 

Throughout this paper, we adopt Q,m ~ 0.3, SIa ~ 0.7, 
(J8 = 0.8, and Ho = 100/ikms"^ Mpc"^ Thus, a.t z = 0.6, 1" 
corresponds to 7AAh~^ kpc and 1' corresponds to 446/i~^ kpc, 
both comoving. At that redshift H{z) = l. 39Ho ,so 5z — 0.1 corre- 
sponds to 216/i~^ Mpc in comoving coordinates. 



2 SAMPLE DEFINITIONS 
2.1 Mg II absorbers 

For the 16713 QSO spectra in the SDSS QSO sample of 
ISchneider svHl j2003h . we searched for Mg II AA2796/2803 ab- 
sorption doublets with Zabs 5S 0.8 using a largely automated tech- 
nique. A third-order polynomial was fitted to overlapping 2500- 
km sections of each QSO spectrum from 10,000 kin s^^ above 
the Lya emission line to 10,000 km s^^ below the Mg II emission 
line. Pixels with flux > 2a below and > 5(7 above the contin- 
uum are rejected and the continuum is re-fitted to the remaining 
points. This process is iterated until no more points are rejected. 
Overlapping portions of adjacent continua are joined by weight- 
ing each linearly from zero at the edge to unity at the centre. The 
final continuum is smoothed over 11 pixels 760 km s"^). The 
2500 km chunk -size is small enough to fit most emission fea- 
tures but large enough so that strong Mg II doublets do not cause 
significant spurious bends in the continuum. 

Candidate Mg II A2796 lines are searched for by identifying 
the pixel within a AArost = 7 A sliding window with the most sig- 
nificant flux deviation > 2 cr below the continuum. This AArcst 
accounts for most of the Mg II A2796 absorption while avoiding 
significant overlap with the Mg II A2803 line. A similar window 
is centred on the Mg II A2803 wavelength with the same redshift, 
2abs, as the putative Mgll A2796 line. If VK.'^sii ^ i.q A and the 
mean signal-to-noise ratio (S/N) is ^ 10 within each window, the 
system is identified as a candidate Mg II absorber. 

Candidates require at least one supporting transition to be con- 
sidered real detections: Fe II A2600 or Mg I A2852. Fe II is the pre- 
ferred line since it is stronger than Mg I. Fe II A2600 is, in principle, 
detectable in SDSS spectra when z^hs > 0.47 (i.e. Aobs > 3800A) 



Table 1. Catalogue of 212 Mgll absorbers from the SDSS DRl with 
0.4 ^ 2abs 5; 0-8. The J2000 name, QSO and absorption redshifts 
and the measured Wr for Mg IIAA2796 & 2803, Mg I A2852 and Fe II 
A2600 ai'e given. Here we show only a small sample from the full ta- 
ble which is available in the electronic edition of this paper and from 
http://www.ast.cam.ac.uk/~mim/pub.html Full name designations and sta- 
tistical uncertainties in Wr are given in the electronic version. 



SDSSJ 



^abs 



Rest equivalent width |^Aj 
2796 2803 2852 2600 



004041-005537 


2.090 


0.6612 


4.37 


3.46 


0.86 


1.83 


004721-1-154652 


1.272 


0.7742 


1.31 


1.92 


0.57 


1.12 


005130-1-004150 


1.190 


0.7394 


2.19 


1.71 


0.58 


1.27 


005408-094638 


2.128 


0.4778 


2.18 


2.09 


1.20 


2.17 



and so we use the criterion advocated by iNestor et alj i2003h to 
select for damped Lya systems (DLAs), ^ O.SA. If Fell 

A2600 is not detectable or if S/N < 10 in the Fell A2600 win- 
dow, we require that Mg I A2852 is detected with S/N ^ 10 and 
W^^^^ ^ O.2A. Several caveats apply to the above requirements, 
particularly when one or more transitions fall near emission fea- 
tures and the broad absorption features often associated with them. 
We will describe these caveats in detail in a later paper, suffice it to 
note here that they apply to < 20 per cent of (real) systems. 

Finally, we remove clearly spurious candidates by visually in- 
specting each Mg II spectrum. The most common mis-identification 
is broad C IV absorption near the C IV emission line. 

With the above algorithm we detected and visually confirmed 
212 Mg II absorbers in the DRl. A typical Mg II absorption doublet 
is shown in Fig.Q Note that, in addition to the crucial lines for 
selection, the Fell A2587 line also often confirms the detection. 
We provide a catalogue of the Mgll absorbers in TableQand the 
absorption redshift distribution is shown in Fig.|2| 



2.2 Luminous red galaxies 

SDSS DRl contains more than 10*^ LRGs over ~ 2000 sq. degrees 
which have lumi nosities Mg < — 21 and fall on the red sequ ence 
with {u-g)o ~ 2 iEisenstein et all200ltlScranton et all2003 V 

For each Mg II absorber, galaxies meeting the following crite- 
ria were extracted from the SDSS DRl galaxy catalogue: 





^pctro 


< 


21, 


(1) 




r*) + 1.2[(r* -r)-0.18] 


> 


1.6, 


(2) 




{g* - r") 


> 


1, 


(3) 


(dx 




> 


0.4, 


(4) 




^psf ^modcl 


> 


0.3, 


(5) 




1 2phot ^atas 1 


< 


0.05. 


(6) 



We also required errors on the model magnitudes to be less than 
0.2 mag in r* and i* , and we excluded objects flagged by SDSS as 
BRIGHT, SATURATED, MAYBE.CR or EDGE. The model mag- 
nitudes were used to compute the colours. Equations are 
the LRG selection criteria of Scranton et al. (2003). Criterion |4| is 
equivalent to imposing Zphot > 0.3. Criterion|5|separates stars from 
galaxies. Criterion |6| is the selection of galaxies within a redshift 
slice of width Wz =0.1 around Z abs using the photometric red- 
shifts, 2phot, of ICsabai et alj JToOSh who showed these to be accu- 
rate to CTz = 0. 1 at r' < 21. The choice of the slice width Wz = 0.1 
corresponds to ~ 200h~^ Mpc (co-moving) and is arbitrary. It is 
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Figure 2. Redshift distribution of tlie LRGs (bottom) and Mg II systems 
(top). The solid curve shows the LRG distribution smoothed with a Gaus- 
sian kernel. The inset shows the distribution of QSO emission redshifts. 



Thus, the amplitude ratio of the cross- to auto-correlation, which is 
('"o.ag/f'o.gg)^ for ^(r) = (r/ro)~^, is a measurement of the bias 
ratio b{Ma_) /b{Mg) which in turn yields the relative halo-masses 
(Afa/Mg). This assumes that ^ag and ^gg have the same slope 7. 

Since our LRG sample is made up of galaxies with photomet- 
ric redshifts, we computed the projected cross- and auto-correlation 
functions, i.e. as a function of physical distance re — 1)^(1 + z)9 
in comoving Mpc with Da the angular diameter distance. From 
the definitions of Wssjre) (e.e.iPhill ipps et al.^ 1 9 7 8,;, Peebles 19931) 
and Wag (rg ) (e.g. lEisensteinllOOSnAdelberger et all2003l) . the am- 
plitude of both u)gg(re) and lUag(ra) is inversely proportional to 
1/Wz, where Wz is the width of redshift distribution For 
a top-hat the ratio uiag(T'9)/ifag(7'6i) is exactly the bias ra- 
tio fe(Afa)/fe(Afg), irrespective of Wz ^. In the case of a Gaus- 
sian redshift distribution Wngire) /wag{rg) is overestimated 
by 25 ± 10 per cent. This factor was determined using (i) numerical 
integration and (ii) mock catalogues (from the GIF2 collaboration, 
lOao et alj|2004l) made of galaxies that had a redshift uncertainty 
equal to the slice width, Wz, as in the case of our LRG sample. 
Note that this factor depends on the shape of not its width. 



a compromise to optimise the signal-to-noise: too small a width 
will yield too few correlated pairs, too large a width will wash out 
the signal. Finally, we remove the 10 per cent of the galaxies with 
problematic photometric redshifts by requiring that galaxies have 
Zphot uncertainties Uz^i^^^ < 0.5. A total of 33,348 galaxies met all 
these criteria in our 212 fields 300 sq. degrees). Fig.|2|shows the 
redshift distribution of these LRGs for the 212 fields. We used the 
spectroscopic redshift when available, which includes only ~ 160 
LRGs. This situation will change in the future with the 2dF/SDSS 
program to obtain spectra of LRGs. 

LRGs are ex pected to have halo-masses > 10^^ M©. 
iBrown et alj i2003h showed that the clustering of red (Bw ~ R> 
1.44) galaxies between z = 0.3 and 0.9 in the NOAO deep wide 
survey is a strong function of luminosity: in the luminosity range 
-21.5 < Mr < -20.5, the correlation length is 6.3±0.5/i"^ Mpc, 
and rapidly increases to 11.2/i^^ Mpc at Mh = —22. Such 
strong clustering is consistent with halo- masses of 3.5 X 10^^ to 
3 X 10" Mq using the bias prescription of lMo & Whit3<2002h . 

3 RESULTS 

3.1 Theoretical background 

A widely used statistic to measure the clustering of galaxies is the 
correlation function, ^(r). The absorber-galaxy cross-correlation, 
^ag, is defined from the conditional probability of finding a galaxy 
in a volume dV at a distance r — |r2 — ri|, given that there is a 
Mg II absorber at ri : 



P(LRG|Mgii) = n„ [1 + eag(r)] dV, 



(7) 



where Uu is the unconditional background galaxy density. 

The observed amplitudes of the auto- and cross-correlation 
functions are related to the dark matter correlation function, ^dm, 
through th e bias, b(M), which is a function of the dark matter halo- 
mass (e.g. lMo et alJl993tlMo & Whitel2002h : 



^gg(r) = fo2(Mg)eDM(r), 
^ag(r-) = 6(A//a)6(Mg)eDM(r). 



(8) 
(9) 



3.2 Mg ii-luminous red galaxy cross-correlation 

Fig- m (filled circles) shows ?j;ag for the entire sample, where 
we used the following estimator of ?j;ag(re) (also advocated by 



we used the rollowing 
lAdelberser et alfcOOsI) : 



/ , AG 
1 + W!,g[re) = -j-^, 



(10) 



where AG is the observed number of absorber-galaxy pairs be- 
tween re — dr/2 and re + dr/2, summed over all the fields. AR is 
the normalized absorber-random galaxy pairs where the normaliza- 
tion is applied to each field independently: AR= AR' /N^, 
where AR' is the number of random pairs in field i, and Ng (N^) is 
the total number of galaxies (random galaxies) for that field. From 
the sample of 33,348 objects within the initial search radius of 40', 
there are 19,496 objects within re = 12.8/i^^ Mpc which is the 
outer radius of the largest bin used. Taking into account the areas 
missing from the SDSS within our search radius, we generated ap- 
proximately 200 times more random galaxies to reduce the shot 
noise of AR to an insignificant proportion. Table |2| shows the to- 
tal number of pairs, AG, and the expected number of pairs, AR, if 
Mg II absorbers and LRGs were not correlated. 

The Wag error bars are computed using the jackknife estimator 
jEfrorl 1982h : we divide the sample into 10 parts and compute the 
covariance matrix from the A'jack = 10 realisations for each part: 



COV» 



jack 



jack 



^ JVjack 

-E 

1 



Wk{rei)-w{rei)]-[wkirej)-w{rej)]{H) 



where Wk is the fcth measurement of the cross-correlation and w is 
the average of the A'jack measurements. 

Two important internal consistency checks on these results 
were performed: using either synthetic Mgll absorbers with real 
LRGs or synthetic LRGs with real Mgll absorbers, we find no 
cross-correlation signal. 



This result is demonstrated in Appendix A. 
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Figure 3. The projected Mgll-LRG cross-correlation, tOag(''e), for 212 
Mgll absorbers and 19,496 LRGs is shown by the filled circles. The er- 
ror bars are computed using the Jackknife technique. The filled triangles 
show the LRG-LRG auto-correlation, Wgg, offset by 0.025 in rg for clar- 
ity, and the dashed line shows the fit to uigg . The solid line shows the rel- 
ative amplitude of tOag and uigg, i.e. using uiag = a X tOgg for scales 
rg > 200/i~^kpc since the smallest scales will be affected by the fi- 
nite cross-section of the absorbers. The best-fitting relative amplitude is 
a = 0.84 ± 0.09. The inset shows the first two bins, 50-100ft,~^ kpc and 
100-200/i~^ kpc (dotted lines), split into two sub-bins. The third sub-bin, 
100-150/i~^ kpc, contains no galaxies due to the absorber cross-section 
effects. 



Table 2. The total number of absorber-LRG pairs AG and the number of 
absorber-random pairs AR, for the cross-correlation shown in Fig. |3] The 
number of pairs expected if one extrapolated u)ag to rg < 0.2h~^kpc are 
shown in parentheses (see text). 



tq AG AR uiag (t(ui) 

[/i-i Mpc] 



0.05-0.1 


4(2.5) 


0.834 


3.797 


2.40 


0.1-0.2 


3(7) 


3.45 


-0.135 


0.625 


0.2-0.4 


20 


14.7 


0.361 


0.274 


0.4-0.8 


72 


57.1 


0.261 


0.146 


0.8-1.6 


263 


232 


0.134 


0.0265 


1.6-3.2 


956 


925 


0.0330 


0.0457 


3.2-6.4 


3702 


3650 


0.0154 


0.0132 


6.4-12.8 


14439 


14200 


0.0171 


0.0281 



3.3 Relative amplitude of cross- and auto-correlation 

In order to constrain the amplitude of lUag with respect to that of the 
LRG-LRG auto-correlation Wgg, we use the following estimator 
with the same galaxies used for the cross-correlation: 

GG 

l + ^«gg(re) = — . (12) 

GG is the total observed number of galaxy-galaxy pairs between 
rg — dr /2 and re + dr /2 and GR is the total galaxy-random galaxy 
pairs, computed as before. The filled triangles in Fig.|3|show uigg. 



The errors and the covariance matrix for w^^ are computed using 
A'jack = 10 jackknife realisations. 

The arguments in section ITTI require that both lUag and Wgg 
have the same slope 7. Therefore, to constrain the amplitude ratio 
Wag/wgg, we first fitted a power law to TOgg(re), and used that as a 
model for Wgg. We will only use the scales larger than 200/i~^ kpc 
in the rest of this paper in order to avoid possible cross-section 
effects (discussed at the end of this section). 

First, the model Wgg{re) = AggVg'^'' for Wgg gives a best 
amplitude, ^gg, at lh~^ Mpc and slope /3gg of 0.175 ± 0.028 and 
—0.780 ± 0.147 respectively ^. The fitted power law Wgg{ro) is 
shown by the dashed line in Fig.|3| 

Then, from the following model for Wag, 

Wag = a X Wgg, (13) 

where a is the amplitude ratio ^ag/j4gg, we find that the best am- 
plitude ratio 

a = 0.84 ±0.09 (14) 

by minimizing « [w — w]"^COV^^[w — w] , where w and 
w are the vector data and model respectively and COV~^ is the 
inverse of the covariance matrix, calculated using single value de- 
composition techniques (see discussion in Bernstein 1994). This 
value is consistent with the fact that 5 of the 6 bins of Wag at 
rg > 200/1 kpc are below Wgg (Fig.|3j- Furthermore, the average 
of the bin ratios Wa,g{re) / Wgg{re) gives 0.85, close to 0.84 ± 0.09. 
Fig.|3|shows both uiag (solid line) and Wgg (dashed line). 

For completeness, a power law fit to w^gire), i.e. Wagirg) = 
Aagrf"*, gives A^g = 0.179 ±0.026, the fitted amplitude at 
Ih-^ Mpc, , and /3ag-0.955 ± 0.118, the slope. 

Note that the relative amplitude, a, is free of systematics from 
contaminants (e.g. stars). This is due to the fact that (1) we use 
the same galaxies for lUag and uigg, and (2) the estimators in equa- 
tionsfTolandfTH are both oc 1 /Ng, where Ng is the number of galax- 
ies. Any contaminants will affect the cross- and auto-correlation 
function in exactly the same way. We find that a is also robust un- 
der numerous different cuts and subsamples. For example, a more 
restrictive star-galaxy separation [equation Js}] gives consistent re- 
sults. Similarly, a is robust to a more stringent cut on the redshift 
difference between the LRGs and Mg II absorbers than in equation 
J6}. Finally, there is no significant variation in a when excluding 
LRGs with larger or smaller SDSS magnitude errors. 

If the Mg II cross-section radius (Rx ) is larger than some of 
the radial bins, this will affect Wa,g{rg) and therewillbea redis- 
tribution of galaxies in the bins near rg ~ fiv.' SteidellJl995 l> con- 
strained Rx to be ^ 70h~^ kpc (comoving) for absorbers with 
W.'^sii ^ 0.3 A. In fact, we find that (1) the first bin at rg = 50- 
100/i~^ kpc is higher than expected if uiag is a single power law 
extrapolated from the large scales (see Table|2}; (2) the second bin 
at rg = 100-200/i~^ kpc is negative and 2 a below the fit in Fig.|3] 
The inset in Fig.|3|focuses on these two bins (indicated by the ver- 
tical dotted lines) which are divided into two smaller sub-bins. The 
region from rg ~ 100/i~^ to 150h~^ kpc (comoving) contains no 
galaxies. Note that since we used projected correlations w{rg), this 

^ The conversion of the amplitude, Agg, to the comoving length, ro, re- 
quires precise redshifts for the LRGs. As mentioned, only ~ 160 specfi'o- 
scopic redshifts are available at present. Nonetheless, a rough estimate is 
ro = 6.2^J^ q/i~^ Mpc using the Mg II redshift distribution and assuming 
a Gaussian ^ with a FWHM Wz = 0.15 for the LR Gs, and is consistent 
with 6.3 ± 0.5/1-1 ^pj, bv lBrownet'Zll2003l) . 
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region corresponds to different angular scales for the range of ab- 
sorber redshifts, 0.4 < Zabs < 0.8. We also find that this signature is 
present even with less restrictive samples. We speculate that given 
the results of Steidel ( 1995), this deficit of galaxies at re — 100- 
150/i~^ kpc is a signature of Rx ■ However, it could be due to some 
other physical process that prevent pairs at that particular scale. In 
a future paper, we will explore these hypotheses with larger data 
sets and simulations. 



4 SUMMARY AND DISCUSSION 

From the SDSS DRl, we selected ~ 20,000 LRGs and 212 Mg II 
absorbers with W^'^^^ ^ lA for Mgll AA2796 & 2803 and 
Vf/f °" > 0.5A for Fe II A2600. We have examined the clustering of 
these LRGs around the Mg II absorbers with unprecedented statis- 
tics on small and large scales (0.05 < re < 13/i~^Mpc). The am- 
plitude of the Mg II-LRG cross-correlation relative to that of the 
LRG-LRG auto-correlation is 0.67 ± 0.07 ± 0.05, after applying 
a correction of 25 ± 10 per cent discussed in section lTTI and in the 
Appendix. The two error terms reflect the statistical and systematic 
uncertainty, respectively. By adding the errors in quadrature. 



a = 0.67 ±0.09. 



(15) 



This corresponds to a correlation length ro,ag 

SMtollh-^ Mpc for Was, and ro 



4.04lQ;7g/i~^ Mpc for the Mg II-Mg II autocorrelation Waa, where 
we used the minimum au to-correlation lengt h ro,gg — 6.3 ± 
0.5h'^ Mpc of LRGs from|Br^net^ 12003). 

Within the context of hierarchical galaxy formation [equa- 
tions {8Jl9)], equation 1151 implies that our Mgll absorbers have 
halo masses 4—20 times smaller than the LRG mininum-mass 
(3.5 X 10^^ Mq), or ~ 2-8 x lO" Mq. For lO" Mq LRG halos, 
the Mg II absorbers have halos of 0.5-2.5 x 10^^ M©. 

Is our mass constrain t of M g II absorbers in agreement with 
the results of lSteidel et alj jT^3) who found that Mg II absorbers 
with W^'^^^ ^ 0.3 A are associated with ~ 0.7L*g galaxies? Our 
mass measurement appears broadly consistent with those results 
given that ~ L*g galaxies have halos of mass ~ IO^^Mq. Further- 
more, the expected amplitude ratio is ~ 0.70, close to our a — 0.67. 
The expected amplitude ratio is found assuming that the correla- 
tion length does not evolve from z — 0.5 to z = and using 
the local correlation o f early and late type galaxies. At 2 = 0, 
IShepherd et alJ j200?) found that the early- and late-type galaxy 
auto-c orrelation lengths wer e ro — 5.45 Mpc and 3.95 respec- 
tively. |Bud^^^^^l]i20o3) found ro = 6.5 Mpc and 4.5 respec- 
tively. Assuming 7 = 1.8 for both of these auto-correlations, then 
from equations [8|c&|9| one expects the late-early c ross-correlation 
amplitude to be (3.95/5.45)^ '*''^ ~ 0.74 (0.72 for lBudavari et alj 
I2OO31) times that of the auto-correlation. 

Note that the re are important diff erences between our Mg II 
sample and that of lSteidel et alj Jl994h . Firstly, our larger equiva- 
lent width threshold, W^^^^ 1 A, will preferentially select sys- 
tems with a larger velocity dispersion over the absorption compo- 
nents. Thus, our sample is potentially biased towards more massive 
halos. Secondly, our M gll sample will be dominated by DLAs: 
iRao & Tumsheld i200Cl) find that ~ 50 per cent of systems with 

^Mgll ^ Q 5 ^ ^FoII ^ Q 5 ^ pj^^g j^ggjpj. 1 2003) 

use W^^^^ 1 A to select a larger proportion of DLAs. 

It should be emphasized that this method (i.e. measuring a 
correlation ratio) has the following advantages: (i) it is free of 
systematics from contaminants (e.g. stars), (ii) it does not require 



knowledge of the true width of the redshift distribution, and (iii) it 
constrains the masses of the Mg II/DLA host-galaxies in a statisti- 
cal manner without directly identifying them. Thus, with a sam- 
ple of confirmed DLAs which are not selected on the basis of 
Mgll line-strength, one should be able to derive the mean mass 
of the DLA host-galaxies with only relatively shallow wide-field 
imaging. This could help establish the relative proportions of low- 
and high-luminosity contributions to DLA host-galaxies. This topic 
is currently under som e debate (e.g. compare iRao et aL.2003i and 
IChen & Lanzettj2003h . 
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APPENDIX A: ON CORRELATION FUNCTIONS^ 

It may seem that taking the ratio between the cross- and auto- 
correlation is inappropriate since the former is based on absorbers 
with spectroscopic (i.e. accurate) redshifts and a sample of galax- 
ies with photometric redshifts (accurate only to ~ 0.1), while 
the latter comprises only galaxies with photometric redshifts. In 
this paper, we have measured the projected correlation function 
'Wp{rg). For a given field (with one absorber) with galaxies dis- 
tributed with ^ , one may think that the auto-correlation is propor- 
tional to J (^)^ dz while the cross-correlation is proportional to 

/ (^) ^ Thus, at first glance, their ratio is therefore not very 
useful. Below we show the situation to be not so trivial. 

First, some definitions and results that will be useful later. For 
a 3D correlation function £(r) = (r/rp) '', th e projected correla- 
tion function ?iip(rp) is tDavis & PeebleJl983h : 



'0 ^7 



(Al) 



where ^(rp, y) is the 3D correlation function decomposed along the 
line of sight y and on the plane of the sky rp, i.e. r^ = j/^ + rp. H^y 
is in fact the Beta function B{a, fe) = t"''^ (1 - t)''^'^ dt eval- 
uated with a = 1/2 and fe = (7 - l)/2, i.e. H-, = ^) = 

r(i)r(V)/r(i). 

In appendix C of lAdelberger et alj i2003h . one finds the ex- 
pected number of neighbours between re ~ dr/2 and re + dr/2 



3 THIS APPENDIX PRESENTS AND ORGANIZES PREVIOUSLY 
PUBLISHED RESULTS, AND IS THEREFORE NOT PART OF THE 
PUBLISHED TEXT. 



6 A^. Bouche, M. T. Murphy, C. Peroux 



within a redshift distance | | < Tz : 



(A2) 



where x = r'^/{r1 + r1) and Ix is the incomplete Beta func- 
tion Bx{a,b) = J^t"-^ (1 - t)'"'^ dt normalized by B{a,b): 
IJa,b) = BJa,b)/B (a,b). 



Many papers jPhillipps et all Il978l : IPeebled Il993t 
iBudavari etM i2003^ have shown that the angular correlation 



function is 



wie) = (9) 



1 — t ^1 



dz{^Yg{z)-\f{zf-- (A3) 



where g{z) — dr/dz — c/H{z) and /(z) = Dc{z) is the comov- 
ing line-of-sight distance to redshift i.e. Dc{z) = dtJ^^. 

Eguation lASI can be derived from the definition s of the angular 
and 3 D correlation functions, w{9) and ^(r) ("e.g. IPhillipps et all 
Il978h . We reproduce the derivation here and extend it to projected 
auto- and cross-correlation functions. The probabilities of finding a 
galaxy in a volume dVi and another in a volume dV2 at a distance 
r = |r2 — ri |, along two lines of sight separated by 9 are 

dP{9) = 7V^ (1^1(1^2 [l + m(6')] or (A4) 
dP{r) = n{zf dVidV2[l + ^{r)] (A5) 

where M is the number of galaxies per solid angle, i.e. dN/ dQ, and 
n{z) is the number density of galaxies, which can be a function 
of redshift. Given that = J n(z)dV{z) and that dV = 
f{z)g{z) dndz,M= jdz^^ jdz n{z)f{z)g{z). 

To relate ■w{9) and ^(r), one needs to integrate equation IASI 
over all possible lines-of-sight separated by 9 (i.e. along z\ and 22) 
and equate it with eauation lA4l 



N'^[l + w{9)] 



d2:i/(2i) g(zi)n{zi) 



dz2f{z2fg{z2)n{z2)[l + C(ri2)](A6) 



In the regime of small angles, the distance r\2 (in comoving Mpc) 
can be approximated by: 



2 

ri2 



+ — 2r\r2 cos 9 
(ri — r2)^ + r^9^ with r ■ 



r-i + r2 



{g{z){zi~Z2)f + f{zf9' with 2 = 

g{zfy^ + f{zfs^ with J/ = 21 - 22 



21 + 22 



(A7) 



wi9) 



Changing variables in eQuation lA6l from (21 , 22) to (2, y), assuming 
the the major contribution is from 21 ~ 22 and using eguation lA?! 
the angular correlation function is 

dzf{zrg{z)Mzr dyg( V/(2)2g2 + g{zYy^) 

[/^°°d2/2(2)<7(2)n(2)]' 

Changing variables io I — g{z)y, using equation lAll and using a 
normalized redshift distribution, i.e. J dz ^ — 1, equation IASI 
becomes 

w{9)= / dz{^Yg{z)-'x{f{z)9t 



■1 J.7 TT 



(A9) 



which leads to equation IA3I (equation 9 in lBudavari et aljl2003h 

and is one version of Limber's equations. 



In this paper, we measured the projected auto-correlation of 
the LRGs, Wgg{ro), where re = f{z)9 Following the same steps 
as above with rg instead of 9, and dV = (drg)^ g{z)dz, Wgg{rg) 



Wgg{rg) = rg 



1 — f ^7 



(AlO) 



In the case of the projected cross-correlation, Wa.g(rg), the 
conditional probability of finding a galaxy in the volume dV2 given 
that there is an absor ber at a known position ri is, by definition 
re.P.lF,isen,stdnl20n3h. 



dP(2|l)(0) = NgdQ2\l + Wp{rg)] 
dP(2|l)(r) = ng{z)dV2[l + i{r)]. 



(All) 
(A12) 



Using the same approximations (equation IA7> and one integral 
along the line of sight 22 (keeping the absorber at 21), one finds 
that the projected cross-correlation is: 



w^g{rg) = 



dz2f{z2)'^g{z2)n{z2)^{ri2) 



(^) aVrl+9{z)Hz,-Z2V) 



dyg{z) {^) p(2)-^C(Vr|+5(2)V) 



(A13) 



1-7 _7 



' 0,ag -'-'7 ) 



(A14) 



where we approximated — with a normalized top-hat of width 
~ , used equation lA2l and the fact that — 1 since a:: ~ 1 
for a typical width Wz of 2QQh^^ Mpc (Section 2.2). Thus, as one 
would have expected, the cross-correlation is inversely proportional 
to the width of the galaxy distribution. Naturally, in equation lAl II 
and IAI2I the redshift of galaxy 1 (i.e. the absorber) is assumed 
to be known with good precision. If the absorber population had 
poorly known redshifts, one would need to add an integral to equa- 
tion lA13l washing out the cross-correlation signal further. This is 
not an issue for our Mg II absorbers. 

For the projected auto-correlation (equation lAlOt . if one ap- 
proximates ^ by a top-hat function of width Wz, then 

/> 00 

^ggire) = {rgf-' rl^^H.x dz g(z) {^f g{z)-^ 



ire) 



1 — 1 't 



H., X 



' 0,gg -'-'7 



Wz X (re) 



1-7 ^7 



' o,gg -'-'7 ) 



(A15) 



which shows that the auto-correlation depends on the redshift dis- 
tribution of the galaxies in the same way as the cross-correlation, 
(Agkoc 1/W^.The reason for this is that the redshift distribution Sfi 
has a very different role with respect to the correlation functions, 
which can be seen by comparing equations lA10l and lA13l It is this 
very different role that leads to the same 1/Wz dependence. 

The above considerations were for one absorber and can be 



■* In general this should be Da (1 + z)9 where Da is the angular distance. 
For a flat universe, Da{1 + z) = Djvf = Dc = /(z) where Dm i s the 
comoving transverse distance, using D. Hogg's notations faogJl99^ . 



The clustering ofLRGs around Mg ll absorbers 7 



easily extended for many absorbers, since the projected correla- 
tions are measured at the same scales (by definition): w-p{rg) = 
° Wp^i{re), where Wp,i is the projected correlation func- 
tion for one field and TVa is the number of absorbers (or fields). 

In the case of a Gaussian redshift distribution the ratio of 
cross- and auto-correlations may not be exactly unity. Using Mock 
galaxy samples (from the GIF2 collaboration, Gao et al. 2004) se- 
lected in a redshift slice of width, Wz , equal to their artificial Gaus- 
sian redshift errors , we find that the cross-correlation is overes- 
timated by 25 ± 10 per cent. Quite importantly, this correction fac- 
tor is independent of the width of the redshift distribution as long 
as — Wz or as long as it is Gaussian. This implies that the ratio 
of the correlation functions {wa.g/wgg) will be insensitive to errors 
in photometric redshifts. 
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